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Abstract. In the last decade there have been many results about special families of graphs 
whose number of perfect matchings is given by perfect or near perfect powers [5] [11] [8]. 
In this paper we present an approach that allows proving them in a unified way. We use 
this approach to prove a conjecture of James Propp stating that the number of tilings of 
the so-called Aztec dungeon regions is a power (or twice a power) of 13. We also prove a 
conjecture of Matt Blum stating that the number of perfect matchings of a certain family of 
subgraphs of the square lattice is a power of 3 or twice a power of 3. In addition we obtain 
multi-parameter generalizations of previously known results, and new multi-parameter exact 
enumeration results. We obtain in particular a simple combinatorial proof of Bo- Yin Yang's 
multivariate generalization of fortresses, a result whose previously known proof was quite 
complicated, amounting to evaluation of the Kasteleyn matrix by explicit row reduction. 
We also include a new multivariate exact enumeration of Aztec diamonds, in the spirit of 
Stanley's multivariate version. 



1. Introduction 

A lattice in the plane divides the plane into elementary regions. A tile is the union of 
two elementary regions that share an edge. A region is a connected region in the plane 
whose boundary consists of lattice segments. A tiling of a region is a way to pair the 
elementary regions it contains into disjoint tiles. 

Tiles are allowed to carry weights, in which case the weight of a tiling is defined to 
be the product of the weights of the constituent tiles. The tiling generating function of a 
region R, denoted T(i?), is the sum of the weights of all its tilings. 

An important motivation for studying tiling generating functions is that tilings of lattice 
regions can be identified with perfect matchings of their dual graphs, and are therefore 
instances of the dimer model of statistical physics on various lattice graphs. 

In the last decade there have been many results about special families of regions on 
various lattices whose number of tilings is given by perfect powers [5] [11] [8]. In this paper 
we present a unified approach that allows proving these cases as well as new results. 

The uniform approach in our proofs is afforded by Theorem 2.3 (the General Comple- 
mentation Theorem) , a result relating the matching generating function of weighted graphs 
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Figure 1.1. The Aztec diamond region of order 4. 



in a large class (that includes in particular all planar bipartite graphs"'^) to the matching 
generating function of their "complement." This result is an extension of Propp's Gen- 
eralized Domino-ShufHing [9], which in turn extended the Aztec diamond case of the 
Complementation Theorem of [2] to arbitrary weights by eliminating a technical assump- 
tion in the statement of the Complementation Theorem of [2] (the assumption that the 
cell- factors are constant along lines of cells — see Section 2 and [2]). Our proofs employ 
combinatorial arguments to bring the problems to a form where Theorem 2.3 applies. 

This more general version can in fact be proved using the arguments employed in [2] 
to prove the Complementation Theorem. However, it is shorter and clearer to prove it 
in the spirit of [9]. We present such a proof, based on certain local graph replacement 
operations, in Section 2. 

A crucial ingredient in our proofs is that certain matrices that encode the type of each 
family of regions behave well under the action of a naturally arising matrix operator d 
whose action is given by certain rational maps (see Section 3). More precisely, in all 
these examples it turns out that after applying d successively to the encoding matrix some 
convenient number of times — in two of our examples, twelve and respectively thirty(!) 
times — , one arrives at a matrix of the same type as the original one, thus obtaining 
simple recurrence relations that determine the tiling generating function. 

In several of our examples we find in fact more general matrices, with entries inde- 
terminates, which also have this periodicity property under d. These can be interpreted 
as parametrized subvarieties of the algebraic variety of matrices having some fixed order 
under the action of d. 

On the square lattice, Elkies, Kuperberg, Larsen and Propp considered in [5] the Aztec 
diamond regions and proved that the region of order n has 2"^"+^)/^ tilings (for n = 4 this 
region is illustrated in Figure 1.1). 

For the square lattice with the diagonals drawn in — the lattice corresponding to the 
afhne Coxeter group B2 — , B-Y. Yang [11] considered regions called fortresses, and proved 
that the number of their tilings is always a power of 5 or twice a power of 5 (Figure 4.1 
illustrates the case n = 4). 

On the triangular lattice, which corresponds to the affine Coxeter group A2 , MacMahon 
[7] showed that the number of tilings of the hexagon of sides a, b, c, a, b, c (in cyclic order) 
is 

" * '' i + j + k — 1 

It is interesting that there is no known family of regions on the triangular lattice whose 
number of perfect matchings is given non-trivially by perfect powers. The "rectangular" 



Indeed, the statement of Theorem 2.3 clearly applies when H is an arbitrary subgraph of the square 
grid. By performing "vertex splittings" (see e.g. [1]) if necessary, any planar bipartite graph can be viewed 
as a subgraph of the square grid, with unaifected luatching generating function. 
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region obtained from a hexagon with sides m, m, n, m, m, n by cutting off the vertices 
not incident to the sides of length n along zig-zag lattice lines has n™ tilings, but it can 
be split into m hexagons that can be tiled independently (see [3]). 

There is only one other afline Coxeter group of order two, G2, and its corresponding 
lattice is the triangular lattice with all altitudes of unit triangles drawn in. Propp con- 
sidered regions on this lattice called Aztec dungeons (see Figure 3.1 for an example), and 
conjectured that the number of their tilings is always a power of 13 or twice a power of 
13. We prove this conjecture and present two weighted generalizations in Section 3. 

B-Y. Yang [11] extended his result on enumerating tilings of fortresses by proving, 
using a quite complicated argument explicit row reduction of the Kasteleyn matrix — , a 
conjecture of Kuperberg and Propp stating that a certain multivariate weighting of the 
Aztec diamond leads to a matching generating function that factors nicely. We present 
a simple combinatorial proof of this general version in Section 4, which also contains a 
weighted form that keeps track of the number of different types of tiles used in a tiling of 
the fortress (for a simple proof of the original version of B-Y. Yang's result, based on the 
Complementation Theorem, see [2]). 

Section 5 contains a couterpart of Stanley's multivariate generalization ([10]; see also 
[2]) of the Aztec diamond theorem of [5]. 

In Section 6 we present a simple proof of an ex-conjecture of Propp, first proved by 
Chris Douglas [4] , stating that the number of tilings of a certain subgraph of the Aztec 
diamond is a power of 2. We also obtain a natural weighted generalization that keeps 
track of the horizontal and vertical edges. 

In Section 7 we consider the tiling of the plane by equilateral triangles, squares and 
regular hexagons. Another ex-conjecture of Propp [8], first proved (but not yet published) 
by Ben Wieland, states that the number of tilings of a certain family of regions on this 
lattice, called dragons (an example is illustrated in Figure 7.1), equals a power of two. We 
present a simple proof of this, and a natural weighted generalization. 

Section 8 presents a proof of a conjecture of Matt Blum stating that the number of 
perfect matchings of a certain family of subgraphs of the square lattice is a power of 3 or 
twice a power of 3. The last section of the paper contains some concluding remarks. 

2. The General Complementation Theorem 

A perfect matching of a graph is a collection of vertex-disjoint edges that form a spanning 
subgraph. For a weighted graph G, the weight of a perfect matching is the product of the 
weights of the constituent edges. The matching generating function M(G) is the sum of 
the weights of all perfect matchings of G. 

We present in this section a generalization of the Complementation Theorem of [2]. We 
review the necessary definitions and notations from [2] below. 

A cellular graph is a finite graph whose edges can be partitioned into 4-cycles — the 
cells of the graph — such that each vertex is contained in at most two cells. 

Let Co be a cell of a graph G and consider two opposite vertices xq and yo of cq. By 
definition, is contained in at most one other cell besides cq. If such a cell ci exists, 
let Xi be its vertex opposite xq. Then xi in turn is contained in at most one cell other than 
ci; if such a cell exists denote by X2 its vertex opposite xi, and continue in this fashion. 
Repeat the procedure starting with yo- 
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The set of cells that arise this way is said to be a line of G. If the sequence {xi} (and 
hence the analogous sequence {yi}) defined above is finite, the line is called a path, and 
the last entry of each of the two sequences is called an extremal vertex of G. 

A weighted graph is a graph equipped with a weight function on its edges. Unless 
stated otherwise, the edge-weights are considered independent indeterminates. Let V{G) 
and X{G) be the sets of vertices and extremal vertices, respectively, of G. Let E{G) he 
the set of edges of G. 

The following definition is a relaxation of the corresponding one from [2] (the original 
definition involved the additional requirement that H is an induced subgraph of G). 

Definition 2.1. Given a weighted graph H, a weighted cellular graph G is said to be a 
cellular completion of H if 

{i) H is a subgraph of G, and all edges of E{G) \ E{H) with both endpoints in H have 
weight in G, and 

(ii) V{G)\V{H)CX{G). 

Therefore, any (not necessarily induced) subgraph of a weighted cellular graph has a 
cellular completion; see Remark 2.4. A graph may have more than one cellular completion: 
Figure 2.1(c) shows a cellular completion of the graph in Figure 2.1(a); another is obtained 
by shading in in Figure 2.1(b) the other class of the chessboard indicated coloring. The 
cellular completion of a subgraph of the square grid is not necessarily a subgraph of the 
square grid. 

Let G be a cellular completion of the graph H. The complement H' of H (with respect 
to G) is defined as follows. First we define the underlying graph of the complement (the 
weight on H' is defined in the next paragraphs): it is the induced subgraph H' of G 
whose vertex set is determined by the equation V{H^)AV{H) = X(G). where the triangle 
denotes symmetric difference of sets (an example is illustrated in Figures 2.1(a)-(d)^; the 
edges of G not contained in H are represented by dotted lines in Figure 2.1(c)). In other 
words, V{H') is the set obtained from V{H) after performing the following operation at 
each end of every path of G: if the corresponding extremal vertex belongs to V{H) , remove 
it; otherwise, include it. 

Starting with the weight wt on G, we define a new weight wti on G as follows (the 
weight of the complement H' will be defined by restriction from wti). If the cell c has 
edges weighted by x, y, z and w (in cyclic order), set A(c) := xz + yw. We call A(c) the 
cell-factor of c. 

Consider a cell c of G. Assume first that the set of endpoints of edges of H contained 
in c is the full set of four vertices of c — we call such cells whole. In this case, for each edge 
e of c set wti(e) := wt(/)/A(c), where / is the edge opposite e in c (this is summarized 
in Figure 2.2, where the left part shows the original weight on a whole cell and the right 
part the new weight). If c is whole it follows that A(c) 7^ 0, and this is well-defined. Note 
that this definition has the effect of sliding each 0-weight edge to the opposite edge in its 
cell, so by interpreting these back as missing edges, the graph changes — see Figures 2.1(a) 
and (d), where 0-wcight edges are indicated as missing. 

If on the other hand the set of endpoints of edges of H contained in c has at most three 
elements — such a cell is called partial — , define wti on c as indicated by Figures 2.3 and 

^In Figure 2.1(d) the graph H' is indeed an induced weighted subgraph, with the missing edges weighted 
by (see the definition of the derived weight wt'). 
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Figure 2.1(a). The graph H. Figure 2.1(b). Whole and partial cells of H. 




•- X • 

Figure 2.1(c). A cellular completion of if. Figure 2.1(d). The graph H'. 





Figure 2.2. Weight change on whole cells. 




xNJW^y l/(2x)\ /l/(2y) 

Figure 2.3. Weight change on partial cells with three vertices. 




Figure 2.4. Weight change on partial cells with two vertices. 



!.4. More precisely, if the aforementioned set has three elements, the partial cell c looks 
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like the left part of Figure 2.3. By our definition of the complement, the missing extremal 
point from cell c will be included in H' , making this a whole cell in H' . Assign it new 
weights as indicated by the right part of Figure 2.3. 

In case the set of endpoints of edges of H contained in c has two elements, define wti 
similarly, using Figure 2.4. 

We note that if a partial cell contains an extremal vertex, this cell doesn't generate a 
whole cell in H'\ in this case just use the corresponding restriction of the right side of 
Figure 2.3 or 2.4. Instances of this occur in Figures 2.1(a) and (d). 

Definition 2.2. Let wt he a weight on the edge set of H. The derived weight wt' on the 
complement H' is the weight obtained by restricting the above defined weight wt'i to H'. 

(The above defined weight is related to, but different from the complementary weight 
defined in [2] and denoted there by the same symbol.) 

The following result is a generalization of Theorem 2.1 of [2]. The main improvements 
are that here the cell-factors are not required to be constant along lines, and that in the 
current set-up H does not need to be an induced subgraph of its cellular completion G 
(see Remark 2.4). 

Theorem 2.3 (General Complementation Theorem). Let H be a weighted graph 
and let G be a cellular completion of H. Let wt be a weight function on the edges of H. 
Then M(i7;wt) can he expressed in terms of the complement H' of H with respect to G, 
weighted by the derived weight wt', as 



where Ci and Cq are the sets of whole and partial cells, respectively. 

Proof. The particular case when H is an Aztec diamond is due to Propp [9] and follows 
by repeated application of Lemma 2.5 below. For the general case we proceed as follows. 
Detach the cells from one another and introduce one new vertex and two new unit- weighted 
edges between any two touching cells; introduce two new vertices and two new unit- 
weighted edges at every vertex of a cell that doesn't belong to any other cell (the latter 
are precisely the extremal vertices of the cellular completion G that are contained also 
in H; for the graph in Figure 2.1(a); this procedure is illustrated in Figures 2.5(a) and 
(b)). By the "vertex splitting" trick (see e.g. [1]), the matching generating function of the 
resulting graph Hi is the same as that of H. 

In the graph Hi there is an opportunity to apply Lemma 2.5 around each whole cell, 
and an opportunity to apply either part (a) or (b) of Lemma 2.6 around every partial cell. 
Let i?2 be the graph resulting from Hi after performing all local replacements prescribed 
by these Lemmas (Figure 2.5(c) illustrates H2 for the example of Figures 2.5(a) and (b); 
H2 is obtained by including the dotted 4-cycles and removing all vertices and edges they 
enclose). 

Corresponding to each extremal vertex of G contained in H wc have a forced edge in 
H2 (these are circled in Figure 2.5(c)). Corresponding to each extremal vertex of G not 
contained in H, we have a new vertex in H2, created by application of Lemma 2.6(a) 
or (b). Therefore, the graph obtained from H2 by removing all forced edges (which by 
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Figure 2.5(a). 




Figure 2.5(b). Figure 2.5(c). 



construction have weight 1) is precisely the complement H' (compare Figures 2.5(c) and 
2.1(d)). By Lemmas 2.5 and 2.6, the resulting weight on H' is precisely wt'. □ 

Remark 2.4. Any subgraph H oi a. weighted cellular graph G has a cellular completion. 

Indeed, restricting G to the cells containing at least one edge of H, redefining as being 
the weight of all edges in E{G)\E{H) having both endpoints in H, and finally "detaching" 
the remaining cells of G at vertices in V{G) \ V{H) (so that they become extremal), we 
obtain a cellular completion of H. Therefore, Theorem 2.3 applies to any subgraph of a 
weighted cellular graph. 

Let G be a weighted graph containing a subgraph isomorphic to the graph K shown 
in the left part of Figure 2.6 (the labels indicate weights; unlabeled edges have weight 1). 
Suppose in addition that the four inner vertices of K have no neighbors outside K. Let 
s(G) be the graph obtained from G (by "surgery") by replacing K by the graph K shown 
on right in Figure 2.6 (dashed lines indicate new edges, weighted as shown). 

The following result is a generalization due to Propp of the "urban renewal" trick first 
observed by Kuperberg (which corresponds to the case when all weights are 1). It follows 
by analyzing the restrictions of matchings of G and s{G) to K and K, respectively (see 
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Figure 2.6. 



yl2/' ~~.,a:/2 



A, . , ,C A< 



V2x\^ / ll2y 
B *B 

Figure 2.7. 

1/2 -'' ;c/2 



l/2;c\_ 1/2 

Figure 2.8. 

[9]). 

Lemma 2.5. M(G) = (a;^ + yw) M(s(G)). 

The next result appears to be new. Together with Lemma 2.5, it is an essential ingre- 
dient in the proof of Theorem 2.3. 

Lemma 2.6. (a) Consider the above local replacement operation when K and K are the 
graphs shown in Figure 2.7, with the indicated weights {in particular, s{G) has one vertex, 
D, that was not a vertex of G, which is incident only to A and C) . We have 



M(G) =2M(s(G)). 

(&) The statement of part (a) is also true when K and K are the graphs indicated in 
Figure 2.8 {in this case s(G) has two vertices, C and D, not belonging to G; they are 
adjacent only to one another and to B and A, respectively) . 

Proof. To prove part (a), partition the perfect matchings of G into three classes: (1) 
those containing the edge of K weighted x and the edge of K incident to C, (2) those 
containing the edge of K weighted y and the edge of K incident to A, and (3) those 
containing the three edges of K incident to A, B or C . Partition the perfect matchings of 
s(G) also in three classes: (1) those containing edge CD but not AB, (2) those containing 
AD but not BC, and (3) those containing two opposite edges of the 4-cycle ABCD (this 
is indeed a partition, since D is incident only to A and C). 
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Figure 3.1. The lattice G2 and the Aztec dungeon D/^. 



Given a perfect matching /z of G, construct a perfect matching /x' of s(G) by discarding 
the edges of [i with both cndpoints in and inchiding edges of K as indicated by the 
correspondence of like-numbered classes in the partitions from the previous paragraph. It 
is easy to check that — > /u' is a bijection and that wt(/i) = 2 wt(/i'). This proves part (a). 

Part (b) is proved similarly. □ 

When H is an Aztec diamond we obtain the following important special case due to 
Propp (called Generalized Domino-Shuffling in [9]). 

Corollary 2.7 (Reduction Theorem [9]). 

M(^£>„; wt) = M(A£)„_i; wt') JJ A(c). 

Provided all cell-factors are nonzero and keep being so after each application of the Re- 
duction Theorem, Corollary 2.7 can be applied successively — with the weight changing at 
every step — until we get down to the Aztec diamond of order 0, whose matching generating 
function is equal to 1. This provides in particular an algorithm for finding the matching 
generating function of the original weighted Aztec diamond. In certain situations, the 
weight pattern repeats after a small number of successive applications. In these cases the 
matching generating function turns out to be given by perfect (or near-perfect) powers. 
Several instances of this phenomenon arc presented in the following sections. 

3. Aztec dungeons 

The plane lattice corresponding to the afflne Coxeter group G2 is obtained from the 
equilateral triangular lattice by drawing in all altitudes (see Figure 3.1). Consider the 
family of all lattice lines parallel to a given lattice line direction and also the family 
of all lattice lines perpendicular to t. The union of these two families forms a sublattice 
isomorphic to the square lattice. 

Tracing out what would correspond to an Aztec diamond region of order n on this 
sublattice we obtain the Aztec dungeon D„ of order n (see Figure 3.1 for an example). 
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Propp conjectured that the number of tihngs of an Aztec dungeon is cither a power of 
13 or twice a power of 13. In this section we prove two different generahzations of this 
conjecture (see Theorems 3.1 and 3.8). 

There are three possible shapes of tiles on the lattice G2'- an equilateral triangle, an 
obtuse triangle and a kite-shaped tile. It is natural to keep track of them by assigning 
them distinct weights. Since the total number of tiles in a tiling is determined by the order 
of the Aztec dungeon, we may assume without loss of generality that one of our weights 
is equal to 1. Assign therefore our three types of tiles weights y, x and 1, respectively. 
Denote by wt the weight obtained this way on Z)„. 

Theorem 3.1. Let P denote the irreducible polynomial P := + 3x^y'^ + 3a;^y^ + 
2x? + 2.xy^ + 1. The tiling generating functions of the first six Aztec dungeons Dn are 
given by 





wt) 


= 1 


Tpi 


wt) 


^x^+y^ 


T{D2 


wt) 


= x^y^P 


T(A3 


wt) 


= x^y^P^ 


T{D4 


wt) 


= x'V\x^ 


T(I)5 


wt) 





and for n> 5 we have the recurrence 

T(Z)„+i; wt) = a;8n-16yl6n-20p4„-8 T{D„_,; wt). (3.1) 

Propp's original conjecture corresponds to the following special case. 

Corollary 3.2. The number of tilings of the Aztec dungeon Dn is given by the equalities 
T(i?o) = 1. T(Di) = 2, T{D2) = 13, T{D3) = 13^, T{Di) = 2 • 13^ T{D5) = 13^ and 

for n> 5 by the recurrence 

T(D„+i) = 13^"-8t(D„_5). 
Proof. Set a; = y = 1 in Theorem 3.1. □ 

The proof of Theorem 3.1 will follow from the following three preliminary results. The 
first of them expresses M(Z?„; wt) in terms of a certain weighted count of perfect matchings 
of the Aztec diamond graph AD2n-2- 

Let ^ be a given k x I matrix with k and I even. The centers of the edges of the Aztec 
diamond graph ADn form a 2n x 2n array. Place a copy of A in the upper left corner of 
this array and fill in the rest of the array periodically with period A (i.e., translate A to 
the right in the array I units at a time and down in the array k units at a time; if 2n is 
not a multiple of or Z some of these translates will fit only partially in the array). 

Definition 3.3. Define the weight wt^ on the edges of ADn by assigning each edge the 
corresponding entry of A in the array described above. 
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Figure 3.2. The Aztec dungeon D4 and its dual graph. 




Figure 3.3. Applying Lemma 2.5 to the graph dual to D4. 
Lemma 3.4. Let N be the matrix 



N : 



r ^ 


y 



X 
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-j,2_j_y2 
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3.2_|_y2 
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1 
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y 


y 












M{AD2n 



(3.2) 



We have 

(3.3) 

Proof. Construct the graph dual to the region i.e. the graph whose vertices are 
the elementary regions of £)„ and whose edges connect precisely those elementary regions 
that share an edge (this is illustrated in Figure 3.2 for n = 4). Weight each edge by the 
weight of the corresponding tile. 

Clearly, the tilings of the region £)„ can be identified with the perfect matchings of 
its dual graph. Therefore, the tiling generating function T(Z)„;wt) is the same as the 
matching generating function of the dual graph of 

This dual graph contains many local configurations like the one described at the begin- 
ning of Section 2, providing one with as many opportunities to apply Lemma 2.5. By their 
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Figure 3.4. A graph embedded in 1?. Figure 3.5. A weight on AD^,. 

geometric orientation, these local configurations can be grouped into three families. The 
largest family contains r? members (these are indicated in Figure 3.3). Apply Lemma 2.5 
at each of these r? places. The resulting weighted graph is readily seen to be isomorphic 
with a weighted subgraph of the square lattice (shown in Figure 3.4 for n = 4). 

Furthermore, because of the vertices of degree one. the resulting subgraph of the square 
lattice will have some edges around the boundary that are forced to be contained in all of 
its perfect matchings (the forced edges are shown in thick lines in Figure 3.4). It is easy 
to see that all these forced edges except 4n of them correspond to kite-shaped tiles, and 
have therefore weight 1. The remaining 4n edges correspond to equilateral triangle tiles, 
thus having weight y. 

Removing all vertices (together with all incident edges) connected by forced edges we 
are left with a weighted spanning subgraph of the Aztec diamond of order 2n — 2 (Figure 
3.5 illustrates the case n = 4; dotted edges have weight 1/2). This can be regarded as the 
entire graph ADin-i by weighting all missing edges by 0. By Lemma 2.5, the obtained 
weight comes out to be exactly wtAr, with the matrix N given by (3.2). Taking into 
account the contribution of the 4n forced edges of weight y and applying r? times Lemma 
2.5 we obtain (3.3). □ 

For a kxl matrix A with k and I even define a new k x I matrix d{A) as follows. Divide 
matrix A into 2x2 blocks 

X w 

y ^ 

and assume xz + yw for all such blocks. Replace each such block by 

z/{xz + yw) y/ (xz + yw) 
w/{xz + yw) x/{xz + yw) 

and denote the resulting k x I matrix by B. Define d(^) to be the k x I matrix obtained 
from B by cyclically shifting its columns one unit up and cyclically shifting the rows of 
the resulting matrix one unit left. 

12 



The next simple observation provides a very convenient way of keeping track of the 
evolution of periodic weights on Aztec diamonds when one applies the Reduction Theorem. 



Lemma 3.5. Let A be a k x I matrix with k and I even and consider the weight wt^ it 
determines on ADn according to Definition 3.3. We have 

wt'^ = wtd(A), 

where the weight on the left hand side is defined as in Definition 2.2. 

Proof. This follows from Definition 2.2 and the above construction of d(A). The rea- 
son we need the cyclic shifts in the construction of d(A) is because the derived weight 
of Definition 2.2 is defined by viewing the Aztec diamond of order n — 1 as embedded 

concentrically into ADn- D 

Due to the location of the O's in the matrix N given by (3.2), all the edges in an Aztec 
diamond that are weighted under wtjv are parallel among themselves. By Definition 2.2, 
this property is preserved by the derived weight. In particular, in any cell at least one pair 
of opposite edges is assigned nonzero weights. Therefore all cell- factors are nonzero and 
the Reduction Theorem can be applied successively. 

By Lemma 3.5, the successive weights that occur are the weights corresponding via 

Definition 3.3 to the iterates d^*^(A^), for i = 1, 2, 3, Therefore, if one of these iterates 

would be the same as N — or indeed the same up to a scalar multiplicative factor — then 
we would get a recurrence for the matching generating function on the right hand side of 
(3.3). By (3.3) this would then translate into a recurrence for M(D„; wt) and would solve 
the problem of computing the latter. 

The computation of the iterates d*-*-* {N) can be done very easily with a computer algebra 
package like Maple. The first few iterates don't look very promising, but perseverance pays 
off: the twelfth(!) iterate turns out to be, up to a scalar multiple, exactly the same as N. 
More precisely, one obtains 

A^^'^\N) = kQN, (3.4) 

with 

y4(x3 + xy^ + l)\x^ + 2x^7/^ + + x)^ 
(x2 + y'^Y{x^ + 3a;4y2 + Zx'^yi + ^6 + 2x^ + 2xy^ + 1)4 ' ^ ' ' 

Since at each application of the Reduction Theorem the order of the resulting Aztec 
diamond decreases one unit, it follows that 

M(A£)2„; wtjv) = ki M{AD2n-i2; wtjv), (3.6) 

where the constant fci is the product of all the cell-factors arising in the twelve applications 
of the Reduction Theorem, multiplied by fc^^""^^)^^" (the latter factor is due to (3.4) 
and to the fact that each perfect matching of ADn contains n{n+l) edges). To carry out 
the computation of fci by hand is a fairly strenuous enterprise, but with the assistance of 
Maple it is quite easy. We obtain by (3.6) the following result. 
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Proposition 3.6. For n>6 we have 

M{AD2n; wW) = a;8"-l<5yl6n-44(^2 ^ y2-)24-12„p4„-8 M(A£)2„-12; wt;v), (3.7) 

where P = + Sx^y^ + Sx'^y^ + y^ + + 2xy^ + 1 . □ 

We are now ready to present the proof of the main result of this section. 
Proof of Theorem 3.1. By Lemma 3.4 we have 

T(Z)„+i; wt) = y^^+\x^ + M(AD2„; wtjv) 

and 

T(D„_5;wt) = +2/2)(n-5)^ M(A£>2„-i2; wt;v). 

Together with (3.7) these two equalities imply (3.1). The stated values for the tiling 

generating functions of the first six Aztec dungeons arc easily checked using (3.3) and 
applying repeatedly the Reduction Theorem to evaluate the weighted perfect matching 
counts of the resulting Aztec diamonds. □ 

Remark 3.7. In fact, it turns out that we can further distinguish between the two 

possible orientations of the equilateral triangle tile by assigning them distinct weights in 
our enumeration. Indeed, keeping the weight y for the up-pointing equilateral triangle tiles 
and assigning weight z to the down-pointing ones, the same calculations as above lead to 
the recurrence 

T(£)„+i; wt) = x^''-^^{yzf''-^^{x^ + Sx^yz + Sx'^y^z^ + y^z^ + 2x^ + 2xyz + 1)^""^ 

xT(D„_5;wt), 

for n > 5 (where wt denotes now the new weight). Similarly, we obtain the values 



T{Do 


wt) 


= 1 


T(£»i 


wt) 


= x'^ +yz 


T(£>2 


wt) 


= x^yzQ 


T{D3 


wt) 


= x^yzfQ'' 


T(Z)4 


wt) 


= x^\yzf{x^ +yz)Q'' 


T(Z?5 


wt) 


= x'\yzY^Q\ 



where Q := x^ -\- dx'^yz + ix^y^z"^ + y^z^ + 2x^ + 2xyz + 1. 

One observes that all above polynomials have the property that the exponents of y and 
z in each of their monomials agree. This implies that the same property holds for each 
monomial in the c;xpansion of T(D„; wt), and therefore each tiling of Z)„ contains the same 
number of up-pointing and down-pointing equilateral triangles. 

However, the ultimate such generalization that would keep track also of the orientations 
of the other two kinds of tiles (six possibilities for each) does not seem to lead to a weighted 
count that is expressible by a simple product formula. Indeed, as the Reduction Theorem 
is applied successively larger and larger irreducible polynomials occur in the factorization 
of the resulting A-factors and their form does not seem to fit a simple pattern. 
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The key fact that ahowcd our proof of Theorem 3.1 is the fact that the weight function 
obtained by successive appUcations of the Reduction Theorem repeated — up to a scalar 
multiple — after a finite number of applications (in our case twelve). 

As mentioned in Remark 3.7, this does not seem to happen if we weight each tile 
orientation by a different weight. In the following result we present a weight function on 
Aztec diamonds that depends on eight free parameters which does have the property that 
the weight resulting at the twelfth iteration is a constant times the original weight. The 
result originally conjectured by Propp corresponds, via the argument that proved (3.3), 
to setting parameters a and b equal to 1/2 and the rest equal to 1. 

Consider the 4x4 matrix 

d d 

9 
h 



M 



and the weight wtM it determines on ADn by Definition 3.3. 

Theorem 3.8. Let R denote the irreducible polynomial R := abgh + acfg + bdeh + 2cdef. 
The matching generating functions of the first twelve Aztec diamonds weighted by wtM are 
given by 



M{ADo 


wtM) 


= 1 


M{ADi 


WtM) 


= de 


M{AD2 


WtM) 


= abR 


M{AD^ 


WtM) 


= {defabghR 


M{ADa 


WtM) 


= 2{abfcdefR^ 


M{AD^ 


WtM) 


= 2{def{abghfcfR^ 


M{ADq 


wtM) 


= 2^{abf{cdeffghR^ 


M(A£)7 


wtM) 


= 2^{def{abghf{cffR^ 


M{ADs, 


wtM) 


= 2^{ab)\cdeff{ghfR^ 


M{ADq 


wtM) 


= 2^{def^{abgh)\cffR^ 


M{ADw 


wtM) 


= 2^{abY^{cdeff{ghfR^^ 


M(A£>ii 


wtM) 


= 2^{de)^^{abgh)^°{cffR^^ 



and for higher orders we have the recurrences 

M(A£>2„; WtM) = 24"-12(a6)4"-10(crfe/)^"-12(5/i)4n-16^4n-8 M(A£>2„-12; 

and 



wt 



M) 



(3.8) 



M(^Z)2n-i;wtM)=24»-16(de)4-10(a65/l)4-14(c/)4n-16^4n-12^(^^^^_^^.^^^)_ 

(3.9) 

Proof. As we noted in the proof of Proposition 3.6, Lemma 3.5 implies that the successive 
weights one obtains when applying the Reduction Theorem repeatedly are the weights 
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corresponding to the iterates d''*''(Af), for i = 1,2,3,.... The twelfth iterate is again a 
constant multiple of M. More precisely, we obtain 

S^'^\M) = k2M, (3.10) 

with 

_ 1 {ag + def{bh + cf)^{ag + 2de)^ibh + 2cf)^ 
^ 16 {abgh + acfg + hdeh + 2cdefY ' ^ ' ' 

Since at each application of the Reduction Theorem the order of the resulting Aztec 
diamond decreases one unit, it follows that M(AZ)„;wtM) differs from M(A_D„_i2; wtM) 
just by a multiplicative constant, equal to the product of all the cell-factors arising in the 
twelve applications of the Reduction Theorem, multiplied by fcg" (since ADn has 

(n — 12) (n — 11) edges in each perfect matching). The explicit value of this multiplicative 
constant depends on the parity of n. (In Proposition 3.6 we were only interested in even 
orders because those are the ones that correspond to Aztec dungeons.) 

With the assistance of Maple the multiplicative constant is easily determined. It comes 
out to be the one shown in (3.8) for n even, and the one in (3.9) for n odd. □ 

Remark 3.9. Recurrence (3.1) explains why there are no new prime factors (in the 
polynomial ring Z[a,y]) entering the factorizations of T(_D„;wt) for n > 6 besides the 
ones occurring for the first six values of n. However, we do not have an explanation for the 
aesthetically pleasing and somewhat mysterious fact that the only prime factors occurring 
for n < 5 are the weight indeterminates themselves and the irreducible polynomial P. 
Indeed, (3.5) illustrates two other irreducible polynomials that occur as divisors of cell- 
factors arising when the Reduction Theorem is applied successively they specialize to 
3 and 5, respectively, for a = y = 1 (and they are the only other prime divisors of 
cell- factors arising in this process). However, as indicated by the initial values stated in 
Theorem 3.1, both these prime factors cancel out in the expressions giving T(D„;wt) for 
n = 0, 1, . . . , 5 and hence the only prime factor of the number of tilings of D„ is 13, 
except for the occurrence of a single factor of 2 for n—\ (mod 3). 

The same remark applies to recurrences (3.8) and (3.9) and the initial values stated in 
Theorem 3.8, as (3.11) indicates four additional potential prime factors — ag + de, bh + cf, 
ag + 2de and bh + 2c/ (again, they are the only other prime divisors of cell-factors in this 
case) — which, however, cancel out in the expressions of M(^Z)„; wtM) for n < 11. 

There is another way to place the Aztec diamond contour on the sublatticc of rectangles 
described at the beginning of this section. Indeed, on can place it so that the resulting 
region has no forced tiles. Denote the region obtained this way by En, we will still call it 
an Aztec dimgeon. Figure 3.6 illustrates the case n — 5. 

Enumerating the tilings of En turns out to be a distinct problem from the case of Dn- 
However, as indicated in the next result, the resulting number of tilings is strikingly similar 
to that corresponding to _D„. The emimeration of the tilings of En follows by a suitable 
specialization of Theorem 3.8, which is thus seen to be a common generalization of the 
enumeration of the two kinds of Aztec dungeons. 

Theorem 3.10. The number of tilings of the Aztec dungeon En is given by the equalities 
T{Eo) = 1, T{Ei) = 2-13, T(^2) = 13^, TiE^) = 13^ T{Ei) = 2 ■ 13^, T(^5) = 13^2, 
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Figure 3.6. The Aztec dungeon E4 and its dual graph. 




Figure 3.7. Applying Lemmas 2.5 and 2.6 to the graph dual to E4. 



and for n> 5 by the recurrence 

t(£;„) = 134"-8t(£;„_6). 

Proof. Consider the dual graph of En (Figure 3.6 illustrates the case n — 4). Apply 
Lemma 2.5 at the (n— 1)^ places indicated by light shading in Figure 3.7, and Lemma 2.6 
at the 4n places indicated in the same figure by a darker shade. By these lemmas, with 
each such application the matching generating function changes by a factor of two. The 
resulting graph is isomorphic with AD2n+i, and the resulting weight is wts where 

-1/2 1/2 1 1- 

1/2 1/2 1 1 

1 10 1' 
.1 110. 

Therefore we obtain that 



T(£„) = 2("+i)' M{AD2n+i;wtB). (3.12) 
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The Reduction Theorem and Lemma 3.5 yield, after combining all cell-factors of the thus 
weighted AD2n+i, that 



M(A£)2„+i;wtB) = 2" -iM(A£)2n;wtd(B)) 



(3.13) 



where 



d{B) 



■ 1 1/2 1/2 

1/2 1 

1/2 1 

. 1 1/2 1/2 



1/2 
1/2 



1 1 



1 



By specializing a = b = g = h= l and c = d = e = f 



1/2, (3.8) becomes 



M(A£)2„;wtd(B)) = 2' 



i60-24n 



134"-8M(A£)2n-12;wtd(s))- 



(3.14) 



Using (3.12)-(3.14), M(i?„) is expressed as a constant times M(^£'2n-i2; wtci(B)). By 
(3.13), the latter can be expressed in terms of M(^£>2n-ii; wts), which is in turn, by 
(3.12), a constant multiple of M(ii'„_6). The multiplicative constant works out to be 
13'*""^ and we obtain the recurrence in the statement of the Theorem. The initial cases 
follow by specialization, using (3.12) and (3.13), from the initial cases of Theorem 3.8. □ 



B-Y. Yang [11] considered regions called fortresses on the lattice obtained from the 
grid lattice by drawing in all diagonals (Figure 4.1 illustrates the fortress of order 4), and 
proved, using the permanent-determinant method [6], that the number of their tilings is 
either a power of 5 or 2 times a power of 5. 

By applying Lemma 2.5 to the dual graph, enumerating the tilings of a fortress amounts 
to finding the matching generating function of an Aztec diamond weighted in a certain 
pattern with weights 1 and 1/2. B-Y. Yang proved in fact a more general result [11], 
included here as Corollary 4.3. 

However, Yang's proof is quite complicated and does not explain why a nice factor- 
ization exists. This makes it inviting to look for a combinatorial proof. A combinatorial 
argument for the weight relevant to fortress tilings can be found in [2] . We present below 
a combinatorial proof of the general case of Yang's theorem. Our proof also explains why 
such a nice factorization is obtained. 

Consider the square matrix 



4. Fortresses 



■ xi yi 1/xi 1/yi xi yi 1/xi 1/yi 
X2 2/2 l/a;2 1/2/2 X2 y2 l/x2 1/2/2 



A = 



X2n 2/2n '^/X2n '^/Uin X2n 2/2n 'i-/X2n l/2/2n 
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The weight wt^ it defines on ADn according to Definition 3.3. is determined by its 
submatrix 

xi yi 

X2 2/2 



-X2n y2r. 

For simplicity of notation, rename the weight wt^ as y-j-. Let Si := {x2i-iy2i-iX2iy2iY^^ , 
i = 1, . . . , n. Define r to be the operator sending T to the (2n — 2) x 2 matrix 



r(T) = 



xi/Si 
XijSi 
X3/S2 
xq/Ss 
X5/S3 



X2n-2/Sn-l 
X2n-z/ Sn-l 
X2n 



Si/yi 
S'2/2/4 
S2/y3 
Ss/ye 
S3/y5 



Sn-l/y2n-2 

Sn~l/y2n-3 
Sn/y2n 



Yang's theorem (together with an explanation for why a nice factorization exists) will 
follow from our next two results. 



(4.1) 



(4.2) 



Theorem 4.1. Forn >1 we have 

n 

M(A£»„;y5.) = M(AD„_i;y,(r)) [] Ai/5i, 

for n even, and 

n 

M(^£)„;y-r) = M(A£)„_i;y,(r)) n^'' 

i=l 

for odd n, where Aj = X2i-iy2i + X2iy2i-i, for i = 1, . . . ,n. 

Using the above result one can easily prove Yang's formula by induction, provided 

the formula has been conjectured explicitly (it is not difficult to conjecture it based on 
computer calculations of the first few cases). However, a conceptual way to obtain the 
formula is by using the following observation. 

Lemma 4.2. For any 2n x 2 matrix T having rows [xi,yi\, i = 1, . . . , 2n, we have 



r(^)(T) = 



X5 

xe 
xr 



y5 
ye 
yr 



X2n-A 2/2TI-4. 
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Proof. This is checked directly from the definition of the operator r (assistance from a 
computer algebra package like Maple makes this verification quite easy). □ 



By Lemma 4.2 we are guaranteed that successive application of Theorem 4.1 will only 
generate a small number of types of cell-factors and quantities Si (all types appear within 
four consecutive applications of Theorem 4.1). We obtain the following result. 

Corollary 4.3. (B-Y. Yang) The matching generating function M(^£)„;yy) can be 

expressed as a simple product in the variables Xi and yi, « = 1, . . . , 2n. 

In our proof of Theorem 4.1 wc will employ the following simple observation. 

Lemma 4.4. Partition the set of edges of ADn into n + 1 classes as follows. View the set 
of edges as forming a 2n x 2n array. For each i — 1, . . . ,n — 1 include the edges in rows 
2i a,nd 2i + 1 in a class. Let the edges in the first row and those in the 2n-th row form 
two more classes. Then every perfect matching of AD^ contains precisely n edges in each 
class. 

Proof. Let /i be a perfect matching of ADn. For each class, consider the set of n vertices 
of ADn incident only to edges in that class. These guarantee there will be at least n edges 
of fj, in each class. The statement follows from the fact that the number of edges in /j. is 
n(n+l). □ 

Proof of Theorem 4-T Suppose n is even. For notational simplicity, we illustrate in 
detail the case n = 4. The arguments generalize with no difficulty to general even n. 

The original weight j/t on the edges of our Aztec diamond AD4, can be written out 
explicitely as wts, where 



-xi/Ai 




.T2?yi?y2/Ai 


?yi/Ai 


xi/Ai 


xxX2y2l Ax 


X4/A2 


X3X4y3/A2 


x3y3yilA2 


2/4/ A2 


X4IA2 


X3X4y3/A2 


X3/A2 


X3X4y4/A2 


a;4j/32/4/A2 


2/3 /A2 


X3/A2 


X3X4:yi/A2 


xe/As 


X5X6y5/A3 


x-^y^ye/ A3 


ya/As 


X(>/A3 


x^xay^l A3 


X5/A3 


XbX&yal A3 


xeysye/ A3 


2/5/ A3 


a^s/As 


XbXeye/As 


-X8/A4, 


x^xsy-r/Ai 


XYyrys/Ai 


ys/Ai 


Xi/Ai 


Xrxsyr/Ai. 



(where as defined in the statement of the theorem, Aj = X2i-iy2i + X2iy2i-i, for i = 
1, . . . ,4). Taking into account the cell- factors, the Reduction Theorem yields 



■xi yi 1/xi 1/yi xi yi 1/xi 

X2 2/2 l/.^;2 1/2/2 X2 2/2 l/-7'2 

2:3 2/3 l/a;3 1/2/3 X3 2/3 I/.T3 

X4, yi I/X4, 1/2/4 X4 yi I/.T4 

X?, 2/5 l/.^;5 1/2/5 Xr, y5 I/.T.5 

x& 2/6 l/a;6 1/2/6 xq ye 1/xe 

X7 2/7 1/2:7 1/2/7 X7 2/7 1/2:7 

.xs 2/8 1/2:8 1/2/8 2:8 2/8 1/2:8 



i/?yil 

1/272 
1/2/3 
1/2/4 
1/2/5 
1/2/6 

1/2/7 

i/ysJ 



Apply the Reduction Theorem. The derived weight is wtc, with 




\X2i-iyii X2iy2i-iJ 
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Group the rows of C into 4 blocks corresponding to the partition of edges of AD^ from 
Lemma 4.4. By Lemma 4.4 we can scale the blocks to obtain 



M(AD3;wtc) = ]jAr-*M(^D3;wtD), 



(4.4) 



where 



'Xi 


XlX2y2 


X2yiy2 


yi 


Xl 


XlX2y2' 


Xi 


xsxm 


xsVsVa 




Xi 


XiXAyz 




xzXiDi 


XiViVi 


2/3 


x^ 


x^Xiyi 


xe 


XbX&Vb 


XbybVa 


2/6 


Xq 


XbX&y^ 




XBXeVe 


xeysye 


2/5 


X5 


XhX&y& 


-Xs 


XyXsVy 


xryrys 


2/8 


Xs 


xrxgyr- 



D = 



This matrix does not yet have the form of B, but can be brought to that form by further 
scaling of the blocks defined above. Indeed, what characterizes a matrix of type B = (bij) 
is that the product hijhij^2 is equal to 1 for all choices of indices for whic;li both entries 
are defined. Matrix D = {dij) does not have this property, but it does have the property 
that the products dijdij+2 are constant — equal to X2i-iy2i-iX2iy2i — within each block. 
Therefore, Lemma 4.4 can be used to bring our matrix to the form of matrix B. With 
Si = {x2i-iy2i-iX2iy2i)^^^ , i = 1, . . . , 4, we obtain by Lemma 4.4 



(4.5) 



i=l 



where 



E 



By (4.3)-(4.5) we obtain that 



xi/Si 


Si/yi 


Si/xi 


yi/Si 


xi/Si 


Si/yi' 


X4/S2 


S2/y4 


S2/X4 


y4,/S2 


X4/S2 


S2/yi 


X3/S2 


S2/yz 


S2/X3 


2/3/<S'2 


X3/S2 


S2/yz 


xe/Ss 


Sa/ye 


S3/X6 


ye/Sa 


X6/S3 


S3 /ye 


X5/S3 


Sa/yn 


S3/X5 


yn/Sa 


X5/S3 


S3/2/5 


X8/S4 


Si/ys 


Si/xs 


y&/Si 


XS/S4 


Si/ys - 



MiADi- wts) = n ^ M(AD3; wtis). 



i=i ' 

For general even n, the same argument proves (4.1). 

The case of odd n is treated similarly. The only difference is that now the exponents in 
the product in (4.3) are {n + l)/2 and (n — l)/2, as opposed to both being n/2 when n is 
even. Relations (4.4) and (4.5) hold with no change for odd n. One obtains the recurrence 
(4.2). □ 

On the square lattice^ with all square diagonals drawn in there are two types of tiles: a 
square tile and a triangular tile. The latter has four different possible orientations. 

As it was the case for Aztec dungeons, if we assign different weights to the five oriented 
tiles the resulting tiling generating function does not seem to be expressible as a simple 
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Figure 4.1. The fortress of order 4. Figure 4.2. The dual graph of a fortress. 
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\fi a 



Figure 4.3. The natural weight of a fortress. 

product of low-dcgrcc polynomials in the weights. However, if the two triangular tiles 
whose hypothenuses face northeast and northwest have the same weight a, and also the 
two remaining types have the same weight 6, the weighted count is a simple product. 

As noted before the statement of Theorem 3.1, there is no loss of generality in assuming 
that the weight of the square tiles is 1. 

Corollary 4.5. The tiling generating junction of a fortress with tiles weighted as de- 
scribed above is a simple product of quadratic polynomials in a and b. 

Proof. The described tiling generating function is clearly the same as the matching 
generating function of the dual graph of the fortress, with weights on edges induced by the 
corresponding tiles (Figures 4.1 and 4.2 show the fortess of order 4 and its dual graph). 
By choosing convenient local graph replacements of the type described in Lemma 2.5, 
computing it amounts to finding the matching generating function of an Aztec diamond 
weighted with the periodic weight indicated in Figure 4.3. By scaling this weight matrix 
by V2 (with the effect of multiplying the matching generating function by a power of 
•\/2), it becomes of the type addressed in Theorem 4.1. Our statement follows then from 
Corollary 4.3. □ 

5. Aztec diamonds with yet another weight 

We have seen in the previous section a periodic weighting of the Aztec diamond that 
factors into small factors. Another such factorization has been found by Stanley ([10]; see 
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also [2]) for the weight pattern 



xi yi xi yi xi yi x\ y\ 

Xl yi X2 2/2 X2 2/2 X2 2/2 



-X2n 2/2n X2n 2/2n X2n 2/2n X2n 2/2n • ■ 

We show in this section that a third weight pattern, namely the one given by the matrix 



N = 



Xl 

X2 



2/1 

2/2 



2/1 

2/2 



Xl 

X2 



Xl 

X2 



2/1 

2/2 



2/1 

2/2 



Xl 

X2 



X2n 2/2n 2/2n •^2™ 2/2ri 2/2n 2;2n 

also leads to a simple product formula. A matrix having the form above is determined by 
its submatrix 

" Xl 2/1 

X2 2/2 



\-X2n 2/2n- 

consisting of its first two columns. For notational simplicity, denote the weight wtAr by 
Let t be the operator sending T to the (2n — 4) x 2 matrix 



t(T) 



1/X4 
1/X3 

1/xe 
I/.T5 



l/a;2n-2 

l/a;2n-3 



Let Ai = X2i-iy2i + X2iy2i-i, for i = 1, . . . ,n 
Theorem 5.1. Forn > 1 we have 



1/2/4 
1/2/3 
1/2/6 
1/2/5 



l/2/2n-2 
l/2/2n-3 



2ra 



(5.1) 



i=i 

i^2.2rt-l 
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for n even, and 



n 2n 

M(^2)„;yy) = 2"-inAi [J x^'^^'^vt''^'^ nADn-2;^,iT)), (5.2) 

^ ^ i7^2,2n — 1 

for odd n. 

Proof. Suppose first that n is even. As for Theorem 4.1, we illustrate the proof for 
n = 4. Our arguments generalize with no difficulty to general even n. The weights on the 
edges of AD^ form the pattern 



B = 



'Xi 


yi 


2/1 


Xi 


Xl 


2/1 


2/1 




X2 


2/2 


2/2 


X2 


X2 


2/2 


2/2 




X?, 


2/3 


2/3 


X3 


X3 


273 


2/3 


X3 


Xi 


2/4 


y4 


Xi 


Xi 


2/4 


2/4 


0:4 


X5 


2/5 


2/5 


Xi 


X5 


2/5 


2/5 


X5 


Xd 


2/6 


2/6 


Xl 


X6 


2/6 


2/6 


xe 


X7 


2/7 


2/7 


Xl 


X7 


2/7 


2/7 


X7 


-X8 


2/8 


2/8 


Xl 


Xs 


2/8 


2/8 


Xs- 



Apply the Reduction Theorem. The derived weight is wtc, with 



C = 



-xi/Ai 


xi/Ai 


2/i/Ai 


yi/Ai 


xi/Ai 


xi/Ai- 


Xi/A2 


Xi/A2 


2/4 /A2 


2/4/A2 


Xi/A2 


Xi/A2 


X3/A2 


X3/A2 


2/3/A2 


y3/A2 


X3/A2 


X3/A2 


X6/A3 


X6/A3 


2/6/A3 


2/6/ A3 


xe/As 


xe/As 




X5/A3 


2/5/A3 


y5/A3 


X5/A3 


X5/A3 


-Xs/Ai 


Xs/Ai 


2/8/A4 


2/8/A4 


Xs/Ai 


Xs/Ai. 



The Reduction Theorem gives 

4 

M(A£)4; wts) = n ^(^Ds; wtc). (5.3) 

i=l 

By Lemma 4.4, we can factor out the denominators along rows to obtain 

4 

M{AD3;wtc) = n M(AD3; wt^), (5.4) 

i=l 

where 



D = 



'Xl 


Xl 


2/1 


2/1 


Xl 


xr 


Xi 


Xi 


y4 


2/4 


X4 


X4 


X3 


X3 


2/3 


2/3 


2:3 


X3 


Xe 


Xe 


2/6 


2/6 


Xe 


Xe 


X5 


X5 


2/5 


2/5 


X5 


X5 


-Xs 


Xs 


2/8 


2/8 


Xs 


Xs- 
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Apply the Reduction Theorem again. All new cell-factors are monomials in the variables. 
We obtain 



(5.5) 



with 



E 



1/(2x4) l/(2y4) l/(2y4) 1/(2x4) 

1/(2x3) l/(2y3) l/(2y3) 1/(2x3) 

1/(2x6) l/(2y6) l/(2y6) 1/(2x6) 

1/(2x5) l/(2y5) l/(2y5) l/(2a;5). 



By Lemma 4.4, the factors of 1/2 can be factored out, with an overall multiplicative 
contribution of 1/2^. By (5.3)-(5.5) we obtain 



M{AD4;zt) = 2''1[A, H x2y,M(AI)2;zt(T)), 



where T is the submatrix of B consisting of its first two columns. For general even n, the 
same arguments prove (5.1). 

The case when n is odd can be treated similarly. The only difference occurs in (5.5), 
where the exponents of Xj and yi are now both (n — l)/2, as opposed to n/2 and n/2 — 1, 
respectively, in the case of even n. Therefore, (5.3) and (5.4) imply (5.2). □ 

Lemma 5.2. For any 2nx 2 matrix T having rows [xi,yi\, i = 1, . . . , 2n, we have 



h(2) 



(T) 



X5 

X6 

X7 



ye 



X2n-4 2/2n-4. 



Proof. This is readily checked by the definition of the operator t. □ 

Corollary 5.3 (Yang [11]). The matching generating function M(j4£>„;zt) can be 
expressed as a simple product in the variables Xi and yi, i = 1, . . . , 2n. 

Proof. By Lemma 5.2 and application of Theorem 5.1 two successive times we obtain 
that, for n > 4, M(^Z)„; zp) is a simple product times M{ADn-4; zf), where T is the ma- 
trix obtained from T by discarding the first four and last four rows. Repeated application 
of this proves the statement of the Corollary. □ 

6. Squares and hexagons 

Consider the lattice of squares and hexagons shown in Figure 6.1, and regard it as an 
infinite graph H. This graph can naturally be embedded in the square grid, as indicated 
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Figure 6.1. 




in Figure 6.2. Draw an Aztec diamond on the grid (see Figure 6.2; dotted lines indicate 

portions of the boundary lying along missing edges). 

A natural question to look at is to find the number of perfect matchings of the portion 
of the embedding of H that is a subgraph of our Aztec diamond (these subgraphs of H 
have a natural definition also on the original squares and hexagons lattice; the example 
shown in Figure 6.2 is isomorphic, after removing forced edges, to the subgraph of H 
indicated in Figure 6.1). Propp conjectured (in an e-mail post to the "domino" mailing 
list in February 1996) that this number is always a power of 2. This was first proved by 
Chris Douglas [4]. 

We present in this section a weighted generalization of this result that keeps track of 
the orientation of the edges in the perfect matching. As in the previous sections, the proof 
will follow by repeated applications of the Reduction Theorem. 

For definiteness, let us draw the Aztec diamond contour so that its westernmost edge 
belongs to a 4-cycle of H lying outside the contour (this is the way the contour was drawn 
in Figure 6.2). Assign weight 1 to the horizontal edges, and weight a to the vertical edges. 
It is easy to see that the weight on ADn arising this way is precisely wt^, with period 
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Theorem 6.1. For alln> 1, M(A£)2n; wt^) = M(A£)2n+i; wt^) = (1 + a2)"("+^). 

In our proof we will employ the following counterpart of Lemma 4.4. 

Lemma 6.2. Partition the set of edges of ADn into n classes as follows. View the set of 
edges as forming a2nx2n array. For each i = 1, . . . , n include the edges in columns 2i — l 
and 2i in a class. Then every perfect matching of ADn contains precisely n + 1 edges in 
each class. 

Proof. Let /x be a perfect matching of ADn- For each class, consider the set of n + 1 
vertices of ADn incident only to edges in that class. These guarantee there will be at least 
n + 1 edges of /x in each class. Since the number of edges in /j. is n{n + 1), the statement 
follows. □ 

Proof of Theorem 6.1. Apply the Reduction Theorem. Since all the cell- factors are 
either 1 or 1 + a^, we obtain 

M(A£)„;wtA) = (l + a2)'=iM(A£>„_i;wtd(A)), (6.1) 

where ki is a positive integer and 

1 a/(l + a2) 1/(1 + a2) a 1 O" 
a 1/(1 + a2) a/(l + a2) 10 1' 

Apply the Reduction Theorem again. All the cell- factors are now equal to 1, so we obtain 



d(A) = 



M(^£)„_i;wtd(A)) = M(^£)„_2;wtd(2)(^)), 



(6.2) 



where 



d(^)(^) 



The cell-factors for wt 



d(2)(A) 



1 a 1/(1 + a2) 1 a/(l + a2) 
a 1 a/(l + a2) 1 1/(1 + 0^) 

are either l + a^ or (1 -|- a^)~^. One more application of the 



Reduction Theorem yields 

M(A£)„_2;wtd(2)(^)) = {l + a'^f^ M(A£)„_3; wtd(3)(^)), 
where k2 is an integer and 

d^^\A) = 



(6.3) 



l/(H-a2) 1 a l + a^ 0/(1+0^) 
a/{l + a^) l + a^ a I 1/(1 -Fa^) 



Lemma 4.4 clearly holds also when partitioning the edges of an Aztec diamond into vertical 
blocks. Apply this version of Lemma 4.4 to the weight pattern wt^a we have on ADn-z- 
Factor out from every third vertical block, starting with the leftmost one, a factor of 
(1 + a^)~^; leave all other blocks unchanged. By Lemma 4.4 we obtain 



M(A£)„_3; wt<j(3)(A)) = (1 + 02)^=^ M{ADn-3\ wts), 
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(6.4) 



Figure 7.1. The Aztec dragon Rq. 



where fcs is an integer and 

„_r l + a^ 1 a l + a^ a(l + a^)" 
~ |_a(l + a2) l + a^ a I 1 + _ ' 

View B as consisting of three 2x2 blocks. By Lemma 6.2, scaling the entries within a 
block has a simple multiplicative effect on the matching generating function. Factoring 
out the constants (1 +a^)~^, 1 and (1 + a^)^^, respectively, the resulting scaled matrix is 
precisely our original matrix A. We obtain by Lemma 6.2 that 

M(A£)„_3; wts) = (1 + a'^f' M(A£)„_3; wt^), (6.5) 

where is an integer. 

It is easy to see that, because of forced edges, M(j4£)2rx+i; wt^) = M(A£)2„;wtA) for 
all n > 1. Therefore, by (6.1)-(6.5) we obtain, after working out the explicit values of 
fci, . . . , fc4, that 

M(^£»2„+i;wtA) = (l + a2)2"M(A£)2„-2;wtA) = (1 + a^)^" M(A£)2„-i; wt^)- 
Repeated application of this proves the Theorem. □ 

7. The Aztec dragon 

Consider the lattice of hexagons, squares and equilateral triangles illustrated in Fig- 
ure 7.1. Propp [8] considered a family of regions i?„ called Aztec dragons on this lattice 
(i?6 is outlined in Figure 7.1) and conjectured that the number of their tilings is 2"^"+^^ 
This was first proved, in work not yet published, by Ben Wieland (as announced in [8]). 
In this section we use the Reduction Theorem to obtain a weighted version of this result. 

There are two types of tiles. Since each hexagon must be paired with a square, and each 
triangle also with a square, the number of tiles of each type is constant across all tilings 
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Figure 7.2. 



The 



region R'q. 




of a given dragon. Thus, our weighted version (Theorem 7.1 below) wiU not concern the 
natural weight in the original set-up, but a certain other weight that is natural once we 
rephrase the problem. 

Let R'^ be the region obtained from i?„ by including n square-triangle tiles along its 
northwestern boundary, as indicated in Figure 7.2. Clearly, these tiles are forced to be part 
of any tiling of the enlarged region. Enumerating the tilings of i?„ is equivalent therefore 
to enumerating tilings of R'^. 

As in the previous sections, consider the dual graph of R'^ (for n = 6 this is pictured 
in Figure 7.3). It is easy to see that this graph can be deformed isomorphically to the one 
illustrated in Figure 7.4. 

Furthermore, by employing the "vertex-splitting" trick (see e.g. Lemma 1.3 of [1]), the 
latter graph is readily seen to have its perfect matchings identified with those of the graph 
in Figure 7.5 (this figure should be regarded as a subgraph of the grid graph; in particular, 
it has vertices at the midpoints of the segments of length 2). 

However, the graph we obtained this way is just the periodic weighting wt^ of the Aztec 
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Figure 7.4. An isomorphic deformation. 



Figure 7.5. A subgraph of ADi2- 



diamond ADor, , with 



A = 



B = 



Therefore, we have 

M{Rn) = M{AD2n;wtA). (7.1) 

One can also regard the obtained graph (shown in Figure 7.5 for n = 6) as a subgraph 
of AD2n- Weight the vertical edges in this subgraph by 1, and the horizontals by a, i.e. 
consider the weight wts on AD2n, where 

"a 1 a 
1 1 

1 a 

1 a 1 

Theorem 7.1. M{AD2n-wtB) = (1 + a2)"("+i). 

Corollary 7.2. T{R„) = 2"("+i). 

Proof. This follows from (7.1) and the case a = 

Proof of Theorem 7.1. By applying the Reduction Theorem twice and collecting cell- 
factors one obtains 

M(A£>2„; wts) = A^"-' M{AD2n-2; wtc), 

where A = 1 + and C = d^^^ (B) is the matrix 

'a 1 a/A 
A 1 
A a 
.1 a 1/A 
30 



1 of Theorem 7.1. □ 



(7.2) 



C ■ 



1/A' 
a 
1 

a/A, 



By Lemma 4.4, we can factor out 1/ A from the first and last rows of C to obtain 

M{AD2n-2\ wtc) = ((l/A)2"-2)" M(A£)2„-2; wtz,), (7.3) 

wliere 



D = 



-aA 


A 


a 


1- 


A 





1 


a 





A 


a 


1 


. A 


oA 


1 


a. 



In turn, by Lemma 6.2 we can factor out A from the first two columns of D to get 

M{AD2n-2;wtD) = (A2"-1)"-' M{AD2n-2; Wtfi). (7.4) 

By (7.2)-(7.4) we obtain the recurrence 

M(A£»2„;wtB) = A2"M(A£>2n-2;wtB). 

Repeated application of this proves the statement of the Theorem. □ 

Remark 7.3. The number 2"^"+^) of perfect matchings of the above subgraph of AD2n is 
approximately equal to the square root of the total number 2^"(^"+^) of perfect matchings 
of the fuU AD2n- 

8. Powers of 3 

We have seen in the previous sections families of (unweighted) graphs whose number 

of perfect matchings are perfect or near-perfect powers of the base 2, 5 or 13. What 
about 3? Matt Blum has considered (as described in an e-mail message to James Propp 
dated December 1997) a family of subgraphs of the square grid for which he noticed (from 
data on concrete cases) that the number of perfect matchings seems to be always a power 
of 3 or twice a power of 3. We prove in this section that this is indeed the case. 

Consider, following Blum, the sublatticc of the square grid lattice showed in Figure 8.1 
and view it as an infinite graph 5*. Draw the boimdary of an Aztec diamond of order n 
on this lattice in such a way that the easternmost edge has an embedded hexagon cast of 
it. Let Bn be the induced subgraph of S spanned by the vertices lying inside or on the 
Aztec diamond boundary (the graph inside the larger contour in Figure 8.1 is -B13). An 
alternative way to view _B„ is as a subgraph of the lattice showed in Figure 8.2, which is 
just a deformation of the lattice considered by Blum (the graph obtained from Bis after 
removing the forced edges is shown in this guise in Figure 8.3). 

Theorem 8.1. Forn > 31 we have 

M(B„) = 3^^" M(B„_3o), (8.1) 

where x^k+i = 4A; - 12, X5k+2 = 4fc - 10 and x^k+s = x^k+A = x^k+b = 4fc - 8, A; > 6. 

The values of M(B„) for n = 1, . . . , 30 are given in Table 8.1. By (8.1) one obtains that 
M{Bn) is either a perfect power of 3 or twice a perfect power of 3, for all n. 
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Figure 8.1. The region B13. 




Figure 8.2. A deformation of the lattice. Figure 8.3. Bis with forced edges removed. 



Proof. The graph i?„ is the periodic weighting of ADn with period indicated by the 
dashed contour in Figure 4.1. Rotating the indicated subgraph by 180° (due to our choice 
of the "base" edge being easternmost), we read off the period to be the 20 x 20 matrix 
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n 




factorization 




71 


T(B„) 


fa<;torization 


1 


1 


1 




16 


531441 


3I2 


2 


2 


2 




17 


43046721 


^16 


3 


6 


2 ■ 3 




18 


10460353203 


321 


4 


6 


2 ■ 3 




19 


10460353203 


321 


5 


6 


2-3 




20 


10460353203 


321 


6 


6 


2 ■ 3 




21 


10460353203 


321 


7 


27 


3» 




22 


7625597484987 


327 


8 


486 


2 ■ 3^ 




23 


5559060566555523 


333 


9 


486 


2-35 




24 


5559060566555523 


333 


10 


486 


2-3^ 




25 


5559060566555523 


333 


11 


486 


2 ■ 35 




26 


5559060566555523 


333 


12 


6561 


3« 




27 


24315330918113857602 


2 ■ 3«' 


13 


531441 


31^ 




28 


79766443076872509863361 


348 


14 


531441 


3I2 




29 


79766443076872509863361 


348 


15 


531441 


3I2 




30 


79766443076872509863361 


348 



Table 8.1. The first thirty values of T(S„). 



A = 



"1 


1 


1 


1 





1 





1 





1 





1 





1 





1 


1 


1 


1 


1" 


1 





1 





1 


1 


1 


1 


1 


1 


1 


1 


1 





1 





1 





1 





1 


1 





1 





1 





1 





1 





1 





1 


1 
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1 


1 


1 


1 





1 





1 





1 





1 
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1 





1 





1 





1 
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1 


1 


1 


1 
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1 


1 


1 





1 





1 





1 





1 





1 








1 





1 





1 





1 





1 


1 


1 


1 


1 


1 


1 


1 


1 
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1 
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1 


1 
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1 


1 





1 





1 





1 





1 





1 
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1 
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1 
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1 
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1 





1 





1 


1 


1 


1 


1 


1 





1 





1 





1 





1 





1 





1 


1 


1 


.1 





1 





1 





1 


1 


1 


1 


1 


1 


1 


1 


1 





1 





1 


0. 



By Lemma 3.5, the periodic weightings arising by successive applications of the Re- 
duction Theorem have periods given by the orbit of A under the operator d defined in 
Section 3. Using a computer algebra program (for instance Maple) one can easily compute 
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d^*^(^) for sTK'ccssivc vahics of i. One finds that all entries of d*-'^^^(^) are either or 
integer powers of 9, with the position of the zeroes matching perfectly the positions of the 
zeros of A. The exponents of 9 form the pattern 



-oo — cxj 



-oo -oo —1 -oo 1 -oo 1 2—1-4 

0-4-12 1 -oo 1 -oo -1 -oo -oo 

-oo -1 -oo 1 -oo 1 2 -1 -4 4 1 

— 1 2 1 — oo 1 — oo —1 — oo — oo 1 — oo 



-1 -4 4 1 -2 -1 -oo -1 -oo 1 -oo - oo -1 - oo 1 

-oo -oo 1 -oo -1 -oo -1-2 1 4 0-4-12 1 - oc 

4 1-2—1 -oo —1 -OO 1 -oo -oo —1 -oo 1 -oo 1 

-oo 1 -oc -1 -oc -1-2 1 4 0-4-12 1 - oo 1 - oc 



(8.2) 

Denote the matrix (8.2) hy B = {bij)i<ij<2o (the position corresponding to a zero of A is 
recorded in B by — oo). 

The form of B suggests the following generalization of our original weight period A. 
Consider the matrix A{q) = {aij{q))i<ij<2o given by 

aij{q)=q'''\ 1 < i, j < 20, 

with the convention that q~°° = 0. By construction. ^4(1) is our original weight period A. 

Applying the operator d successively thirty times to A{q) (a simple task with the aid 
of Maple), we obtain 

d(3°)(^(g))=A(99). (8.3) 

Therefore, by the Reduction Theorem and Lemma 3.5 we obtain that M(^Z)„; wt^(q)) is 
equal to M(A£)„_3o; wt^jg^)) multiplied by a product of cell-factors. The latter product 
is obtained multiplying together all cell- factors from the (thirtyfold!) application of the 
Reduction Theorem. With the assistance of Maple it is not hard to find this product 
explicitly. One obtains 

M(^£>„; wt^(q)) = 3^^" M(A£>„_3o; wt^(9<,)), n > 31, (8.4) 

where ywk+i = Sk - 13, ywk+2 = yiofe+6 = 8fc - 9, yiofe+3 = ywk+7 = 8fc - 7, t/iofe+4 = 
yiok+5 =8k-8, yiok+s = 8fc - 5 and ywk+g = ywk+io = 8fc - 4, for fc > 3. 
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m(/lUn,^^A(q)) 


fcictoriz3,tioii 








f cictrOr i z 9.tr ion 


1 




9^ 




16 


531441 g2 


3I2 g2 


2 


2g-2 


2g-2 




17 


43046721 g2 


3I6 g2 


3 


6q-2 


2 -39-2 




18 


10460353203 


321 q2 


4 


6 


2 • 3 




19 


10460353203 


321 


5 


6 


2 ■ 3 




20 


10460353203 


321 


6 


6 g 


2 • 3 ij'^ 




21 


10460353203 q^ 


o21 2 

3 g 


7 


27 g2 


33 g2 




22 


7625597484987 q-2 


327 g-2 


8 


486 


2 • 35 g2 




23 


5559060566555523 g-2 


333 q-2 


9 


486 


2 ■ 35 




24 


5559060566555523 


333 


10 


486 


2 ■ 3^ 




25 


5559060566555523 


333 


11 


486 


2 ■ 35 g2 




26 


5559060566555523 q^ 


333 g2 


12 


6561 


38 9-2 




27 


24315330918113857602 q^ 


2 • 340 92 


13 


531441 g-2 


3I2 g-2 




28 


79766443076872509863361 q^ 


3« g2 


14 


531441 


3I2 




29 


79766443076872509863361 


348 


15 


531441 


312 




30 


79766443076872509863361 


348 



Table 8.2. The first thirty values of M(AZ)„; wt^(,)). 

Recurrence (8.4) determines M(A£'„; wt^(q)) for all n provided we know it for n < 30. 
This can be easily found using Maple to carry out successive applications of the Reduction 
Theorem. One obtains the values in Tabic 8.2. These clearly specialize to the values in 
Table 8.1 when q = 1. It is straightforward to check that for g = 1 recurrence (8.4) and 
these initial values imply (8.1). □ 

9. Concluding remarks 

As mentioned in the proof of Theorem 6.1, a result analogous to Lemma 4.4 holds for 
the columns of the edge- weight array. It is easy to see that together with Lemma 6.2 this 
analog of Lemma 4.4 implies that single columns of the cdgc-wcight array can be scaled 
with a simple multiplicative effect on the matching generating function of a periodically 
weighted Aztec diamond. The same is true for the rows of the edge- weight array. 

A natural equivalence relation arises this way on the set of period matrices. In order 
to obtain a recurrence for the values M(AZ)„; wt^), it is enough to find an iterate d^''\A) 
in the equivalence class of A. 

Sometimes it is useful to have an alternative approach, as illustrated in the proof of 
Theorem 8.1, where it was crucial for the proof that we first generalized the original 
weight — by introduction of the parameter q. This larger class of matrices had the desired 
property that a convenient iterate (the thirtieth) was in the same class, yielding thus a 
recurrence. 

All the problems treated in this paper led, in their unweighted form, to periodic weight- 
ings of Aztec diamonds with periods given by certain 0-1 even by even matrices (except 
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for the case of Aztec dungeons and fortresses, when entries equal to 1/2 also arise). A 
natural problem to consider would be to classify the 0-1 matrices that give rise to nice 
enumeration results when regarded from the point of view of the approach described in this 
paper. Two ways for this to come about are mentioned in the previous two paragraphs. 
Given the similarity of the examples considered in Sections 6 and 8, one natural starting 
point would concern sublattices of the square grid obtained by deleting parallel edges in 
some more general simple pattern. 
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